Coexistence of scale invariant states in a scale free system 
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Soft solids, such as elastomers, have an instability toward the nucleation and growth of sharply 
creased surface folds. We develop an analogy between this instability and a first order phase tran- 
sition by constructing an isolated fold as a domain wall-like structure separating two, coexisting 
scale invariant deformations of an incompressible elastomer: an affine deformation, and a crease. 
We relate the existence of this surface folding instability and the striking non-existence of a phase 
boundary separating the two scale invariant deformations to certain convexity properties of the 
elastomer free energy. 



Liquefaction of argon, nematic ordering of the liquid 
crystal 5CB, and allotropy in iron are well studied exam- 
ples of first order phase transitions in materials. These 
phenomena can be understood at a macroscopic scale 
as coexistences between different scale free, homogenous 
states — phases. Recently, soft solids such as a hydro- 
gels, elastomers, and tissues, were discovered to possess 
a novel kind of scale free instability with an uncanny, 
though imperfect resemblance to a first order phase tran- 
sition [IH3] • When a soft solid surface is sufficiently com- 
pressed, infinitesimal, sharply creased folds can nucleate 
and grow [TJ |3j . This critical compression is independent 
of the sample shape and marks the coexistence of local- 
ized folds, called sulci, and a smooth free surface. With 
further quasi-static compression, work flows in from the 
displaced boundaries of the sample to the individual sulci 
as these grow. During growth, the compression close to 
each sulcus remains at the coexistence value [2. The 
amount of work is proportional to the folded volume, and 
hence can be interpreted as the work of transformation. 
Beyond a second, larger critical compression the flat free 
surface loses metastability and sulci spontaneously form 
[TJ 0]. This second critical compression can be under- 
stood as a kind of spinodal point for the flat surface 
[2] The analogy with a phase transition fails, how- 
ever, when we attempt identify a sulcus with a "droplet" 
of a new phase, since, for example, no clear phase bound- 
ary divides the sulcus from the surrounding elastomer 
[2J. In this Letter we propose an explanation for the 
partial analogy between the sulcification (or creasing) in- 
stability of an elastomer surface and a phase transition. 
Whereas phase transitions are transformations between 
distinct homogenous, scale invariant states, instabilities 
like sulcification occur when a system transitions between 
scale invariant, though potentially inhomogeneous states. 
We explain how to construct these states, the transitions 
between them, and the striking absence of phase bound- 
aries. 

Recent interest in sulcification has focused on the ef- 
fects of swelling [BH5], growth [TtM2] . applied fields [15] . 
mechanical confinement [1, 14-18_, and imperfections [19] . 
on how patterns form[2J |3j 120] . and on connections to 
other phenomena such as plastic folding [21 . In previous 



work we calculated the critical strain for sulcificiation [1 
and developed a bifurcation and pattern formation the- 
ory based on treating sulcification point as novel kind of 
nonlinear critical point [2J. The later work was built on 
the supposition that an isolated sulcus in a critically com- 
pressed rubber half-space is a local minimum of energy up 
to translations and changes of scale. At the critical point, 
the scaling mode changes stability, allowing a nucleated 
sulcus to grow from infinitesimal size. Linearization fails 
to detect this instability because scaling changes the size 
of the sulcus without changing the magnitude of its de- 
formation gradient. 

The possibility of such an undetectable small ampli- 
tude instability involving a continuous deformation was 
actually first anticipated by Weierstrass when he distin- 
guished so-called weak (i.e. merely linearly stable) lo- 
cal minima from strong local minima. The detection of 
Weierstrass instabilities has remained an open problem, 
however [22] . Our exposition here adds to these ongo- 
ing discussions by showing how to construct an isolated 
sulcus solution in a simple elastomer, and how this con- 
struction might generalize to more complex systems. 

Our plan is to systematically develop the analogy be- 
tween sulcification and a phase transition. To frame this 
discussion, we consider a sample of incompressible neo- 
Hookean elastomer. This is a model material for many 
rubbery solids which are dramatically softer in simple 
shear than in volumetric compression. It is compactly 
specified by the Lagrangian density 

£(A) = ^A?-p(detL4L]-l) (1) 

where the matrix A = <9x/9X is the gradient of the de- 
formation x(X) of a reference volume. The first term in 
([!} governs the entropic elasticity of the network chains 
with shear modulus fi, and the second term, involving 
the Lagrange multiplier p, enforces the nonlinear incom- 
pressibility constraint det 9x / <9X = 1 . When restricted 
to isochoric deformations so that the constraint term in 
([T| is identically zero, C = T , the free energy density of 
the elastomer. In terms of T , we impose the volumetric 
constraint by setting .F(A) = +oo for dct(A) ^ 1. 

The Lagrangian density C is invariant under rigid 
translations and rotations of x and X, separately, i.e. 
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x — > Ox + c or X — > OX + c where O is an orthogonal 
matrix and c is a vector. The first of these symmetries 
is the so-called objectivity of the model while the sec- 
ond expresses that the elastomer is a homogenous and 
isotropic material. Our focus is on another symmetry of 
C. Because the Lagrangian density £ is a function of the 
dimensionless deformation gradient A, it is also invariant 
under the scale transformation {x, X} — > {Lx, iX} for 
any scaling factor L > 0. Unlike rigid translations and ro- 
tations, however, the measure dX is not invariant under 
scale transformations. Thus while a scale transformation 
is not a symmetry of the Lagrangian L — J n CdX. — and 
so not associated with a Noetherian conservation law — it 
is nevertheless a symmetry of the Euler-Lagrange equa- 
tion. 

Many material models possess scale symmetry as an 
intermediate asymptotic property between the molecular 
length scale and the system scale. The configurations of 
small enough patches of material, which are still larger 
than the small-scale cut-off, must be independent of both 
the macroscopic and microscopic details. In particular 
such configurations must become invariant under rescal- 
ing as the system size is taken to infinity. 

To compare the elastomer with other materials, it is 
helpful to view the deformation gradient A as a tensor 
order parameter. Then, what distinguishes the elastomer 
from other materials with tensor order parameters, say 
nematic liquid crystals, is the compatibility constraint 
that A = 9x/<9X. A similar constraint holds for the 
scalar phase variables that appear in theories of smectic 
liquid crystals and superconductors. What distinguishes 
the elastomer from smectics and superconductors is the 
vector nature of the "phase" variable x. Because of the 
compatibility constraint, an external field conjugate to A 
is equivalent to the application of tractions to the bound- 
ary of a sample. The compatibility constraint plays a 
fundamental role in our analysis below. 

Even though rubbers are stable enough materials to 
use in car tires, let us for the moment forget this stability 
and compute the "phase diagram" corresponding to (JT|), 
viewed as an abstract material. For simplicity, we con- 
sider only plane strain deformations so that x, X £ M 2 . 
In this exercise, we would identify different phases with 
homogenous values of A corresponding to distinct local 
minima of the free energy density T . This is, for exam- 
ple, how the austenite and martensite phases of shape- 
memory alloys are described. If T has multiple local 
minima, then we can consider the coexistence of any two 
phases. Because of the compatibility constraint, phase 
boundaries between such phases can involve at most a 
rank-one jump in A. That is if Ai and A2 are the two 
phases then Ai — A2 = a ® N where N is the inter- 
face normal and a is some vector. (Otherwise the elas- 
tomer material would tear along the phase boundary.) 
One can check that T is rank-one convex|28 . Invoking 
Maxwell's famous argument, we conclude that no two 



phases can coexist for any value of applied field (i.e. im- 
posed stress). Hence the neo-Hookean material has no 
phase transitions. 

Since Maxwell's argument is based on convexity, it is 
perhaps troubling that the neo-Hookean free energy den- 
sity is not actually convex in A. This is easy to see if one 
imagines rotating a homogeneous volume by 180° so that 
A — > — A. The chord joining these states passes though 
A = 0, which is a point of infinite energy because of the 
isochoric constraint. To see the consequence of this non- 
convexity of C, we reframe our exercise to examine an 
infinite elastomer- vacuum interface. We are now looking 
for interfacial phase transitions. 

First let us fix a Cartesian coordinate systems 
(X , X 2 ) in the material frame and (a; 1 , x 2 ) in the labo- 
ratory frame so that the elastomer occupies the material 
space X 2 < 0. We can freely assume the deformation 
x(0) = so that the deformed elastomer also occupies 
and laboratory space x 2 < 0. Because of the free bound- 
ary condition on X 2 = we must have A = diag(A, 1/A), 
where 1 — A is the lateral compression of the elastomer 
interface. We can check the stability of the interface by 
computing the speeds of shear and Rayleigh waves, which 
must be real in a stable material. The shear wave speed 
is always real (because of rank-one convexity), but the 
Rayleigh wave speed is imaginary for compressions ex- 
ceeding 1 — Xb ~ 0.45 [I]. This surface instability is 
known as Biot's instability. Note that if the free energy 
density J- were convex, Biot's instability would not exist. 

In many materials, the vanishing of a wave speed indi- 
cates a phase transition. For instance, the speed of sound 
vanishes at the liquid- vapor critical point. In our present 
model, the shear wave speed vanishes as fi — > 0, and the 
/j = limit of Lagrangian density describes a perfect, 
incompressible fluid. 

Let's, then, conduct the following thought experiment. 
Suppose a thermal fluctuation (or other perturbation) 
triggers Biot's instability. Recent post-buckling analy- 
sis of Biot's instability has shown that Biot's instabil- 
ity is sub-critical, which suggests the hypothetical sur- 
face phase transition is first order [23]. Like most first 
order phase transitions, we expect that the fluctuation 
nucleates an initially localized, growing "bubble" of some 
highly deformed configuration. Because the neo-Hookean 
material has no bulk phase transitions, we might guess 
that the emerging phase in this bubble involves some kind 
of frozen surface wave or roughness. If the amplitude of 
this roughness were finite, we could use the scale symme- 
try of C to rescale the amplitude to zero and recover the 
initial, flat interface at large scales. But as the flat in- 
terface is unstable, it must be that the alternative phase 
has infinite amplitude fluctuations. 

The signature of such a fluctuation-dominated state 
in a finite size system would be the non-attainment of 
the energy minimum. That is, a minimizing sequence of 
deformations starting from a super critically compressed 
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FIG. 1: A sulcus in a unit-diamater half-disc compressed 
to a width A (blowups on right) is mapped to a domain wall 
in a strip (representative configurations A — C) by conformal 

change of variable and coordinates [double headed arrow, 
also Eqs. |3j]. Color (blue to red) indicates a linear blow-up 
of pressure as s — > — oo. The size of the sulcus L is dual to 
the position of the domain wall on the s-axis in the strip. 
The bifurcation diagram (blue line) relates L 2 oc A — A* as 
was explained in [2]; the critical compression 1 — A* as 0.35 

flat interface would develop wilder and wilder surface un- 
dulations at finer and finer scales. A precise criterion 
for the non-attainment of an energy minimum was given 
by C. B. Morrey who proved that so-called quasicon- 
vex energy densities always attain their energy minimum 
[24J|29]. The free energy density T satisfies a stronger 
condition called polyconvexity [2"5| |30| , so its form rules 
out any roughening of the interface. 



Because of the quasiconvex form of J 7 , at long times, 
any finite volume of elastomer surrounding the initial dis- 
turbance must settle into a new equilibrium, which in 
an infinite system must be scale invariant. So we now 
propose to identify all deformations which are invariant 
under scale transformation. 

Scale invariant deformations are most easily identified 
by picking an origin and changing to logarithmic-polar 
coordinates (which is a conformal coordinate transfor- 
mation from the standard Cartesian coordinates on M 2 ). 
We define new coordinates systems (T, 9) and (s, <j>) by 
the relations 

X 1 = cos(e)e T X 2 = sin(9)e T (2) 

x 1 = cos(<j))e s x 2 = sm(</>)e s . (3) 

(I.e. T = logi?, the radial coordinate in the reference 
body, and likewise for s.) This coordinate transforma- 
tion maps an annulus in the the Cartesian coordinates to 
a finite strip in the logarithmic-polar coordinates as illus- 
trated in Fig. [l] In particular, the X 2 < half-plane is 
mapped to the infinite strip (T, 8) E [-co, oo] x [—§,§]. 
More importantly, this coordinate transformation maps 
rescaling in the Cartesian coordinates to simultaneous 
translation in s and T. The coordinate transformation 
also maps power-law type singularities in the deformation 
gradient to T-translation invariant deformation gradients 
on the infinite strip. 

In the original Cartesian coordinates, the Euler- 
Lagrange equation for is dSf/dX a = 0, where the 
nominal stress S" — dC/dA l a . The divergence form of 
this equations is direct consequence of the compatibility 
constraint on A. In the new coordinates, the equations 
of equilibrium take the form (in units where fi = 1) 
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the isochoric constraint becomes 
ds d<j> ds d(f) 
dfdQ ~ dQdf ~ 
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and boundary conditions on the surfaces 9 = ±? (i.e. 
the erstwhile free surface X 2 = 0) are 
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Scale invariant deformations, that is T-homogenous so- 
lutions to Eqs. (5][7l, can be parametrized as ds/dQ = 



f(Q), d(j>/dQ = g(Q), ds/dT = 1, and d<j)/dT = 0. In- 
spection of Eqs. (|5][Tj) shows that such a form requires 
p = CT + h(Q). Solutions with C = correspond 
to afhne deformations. In the logarithmic-polar coor- 
dinates these are s = T + \ log (A 2 sin 2 6 + ^2 cos 2 6) , 
A 2 tan<^> = tan 9. Solutions with C ^ correspond to 
creased deformations with different opening angles j3 and 
with different amounts of shear applied to the bound- 
aries 9 = ±7r/2. All of these deformations develop infi- 
nite compressive stresses normal to these boundaries as 
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T — >• — oo, which would violate the second of Eqs. ^ 
unless /3 = so that the crease closes to self-contact. 
The only solution with (3 = and no surface shear is 
s = T - | log 2, <f> = 26, C = §, and ft = const. |2B]. 
Note that the creased deformation spontaneously breaks 
the ad-translation symmetry of the interface. 

The creased (C ^ 0) deformation and afhne (C = 0) 
deformations are all linearly stable, but for different rea- 
sons. The affine deformations are linearly stable (for 
A > Xb) because the local curvature of the free energy 
density function provides a restoring force for small de- 
formations. The crease is stable because it generates its 
own confining pressure. Formally, this existence of this 
pressure can be traced the terms in Eqs. (|5| which are 
quadratic in s and <f> and to the incompressibility con- 
straint. These terms arise in the change of variable to 
logarithmic-polar coordinates because of the tensor char- 
acter of the nominal stress Sf = dC/dA z a , and the di- 
vergence form of the Euler-Lagrange equation dSf / dX a . 
Thus the same features of the elastomer that yielded the 
non-convexity of T (incompressibility) and the stability 
of the homogenous phase in spite of this non-convexity 
(the compatibility constraint) also stabilize the inhomo- 
geneous scale invariant deformation (i.e. the crease). The 
crease and the affine deformations are distinct manifolds 
of linearly stable, T-translation invariant configurations. 
So just as for two distinct stable phases, there must be 
a domain wall type solution which interpolates from the 
creased manifold at T = — oo to the affine manifold at 
T = +oo. 

Coexistence between the two scale invariant states oc- 
curs when the domain wall can freely traverse the strip, 
see Fig. [I] Transforming this mobile domain wall back 
to Cartesian coordinates yields the sulcus deformation 
of the critically compressed half-space. The location of 
the domain wall in the strip gives the size of the sulcus 
in the half-space, hence the mobility of the domain wall 
is equivalent to the soft growth mode of a sulcus in a 
critically strained half-space. While the width of most 
domain walls is set by microscopic physics, the rank-one 
convexity of T precludes the formation of sharp phase 
boundary. Hence the width of our domain wall must be 
fixed by the width of the strip. The domain wall turns 
out to have an exponential profile, which translates to 
power law tail of the sulcus in the Cartesian coordinates 
and gives rise to the energy of transformation |2J . 

Practically, we find the domain wall solution by solving 
Eqs. (5j7) numerically using finite element and continu- 
ation methods as in Ref. pQ . We impose the creased and 
affine deformations as boundary conditions on opposite 
ends of a long strip in the T — 6 plane. The two imposed 
deformations can have a relative shift in the x 2 direction. 
The value of this shift, which is the depth of the crease 
within the sulcus, is fixed (relative to the size of the sul- 
cus) by requiring that the integral of the traction in the 
x 2 direction vanishes on either constrained surface. We 



find coexistence between the crease and an affine defor- 
mation at the critical compression of 1 — A* m 0.35, which 
is similar to what has been reported elsewhere [U [15] . 

The principle implied our discussion of the incompress- 
ible neo-Hookean model is that while ordinary phase sta- 
bility follows from rank-one convexity of the free energy, 
this property is not strong enough to preclude all pro- 
cesses which resemble phase transitions. In general, scale 
free material models can have multiple scale invariant po- 
tential wells, and some of these may correspond to inho- 
mogeneous states. However, a local analysis of the free 
energy density will only detect the homogenous states, 
i.e. phases. Rather, one should look for all scale in- 
variant configurations, whether or not these are homoge- 
nous. Whenever a model has multiple scale invariant 
potential wells, whether or not they correspond to dis- 
tinct phases, there are associated coexistence conditions 
between every pair. When coexistence occurs, the mode 
in which one state transforms into the other (i.e. the 
growth mode) is an emergent soft mode of the whole sys- 
tem which can dramatically change its response to exter- 
nal stimuli. It is interesting to note that our construction 
doesn't directly use the free energy density J 7 , so it makes 
sense even for systems which are out of equilibrium. We 
also remark that integrability of the strain energy den- 
sity close to the origin X = implies the coordinate- 
transformed deformation gradient must in any case be 
bounded as T — > — oo. So a potential alternative to a 
T-translation homogenous state, such as a crease, is a T- 
periodic state. Such configurations might be relevant in 
the finely twinned microstructure that sometimes forms 
during a martensitic phase transition. Our method of 
constructing a sulcus potentially has an immediate appli- 
cation in liquid crystal elastomers. A recent experiment 
suggests that sulcification or creasing in these materials 
may induce topological changes in the nematic director 
field [27]. The existence of a director field in these ma- 
terials potentially enriches the space of inhomogeneous, 
scale invariant deformations beyond the crease studied 
here. 
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